Are there really phase transitions in 1-d heat conduction models? 
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Recently, it has been claimed (O. V. Gendelman and A. V. Savin, Phys. Rev. Lett. 84, 2381 
(2000); A.V.Savin and O.V.Gendelman, arXiv: cond-mat/0204631 (2002)) that two nonlinear clas- 
sical 1-d lattice models show transitions, at finite temperatures, where the heat conduction changes 
from being finite to being infinite. These are the well known Frenkel-Kontorova (FK) model and 
a model for coupled rotators. For the FK model we give strong theoretical arguments why such 
a phase transition is not to be expected. For both models we show numerically that the efi'ects 
observed by Gendelman et al. are not true phase transitions but are rather the expected cross-overs 
associated to the conductivity divergence as T ^ and (for the FK model) T — > cxd. 



Heat conduction in classical one-dimensional lattices 
has recently been investigated by many authors (see 
for a recent review). A large class of 1-d systems can be 
described by the general Hamiltonian 



/ 2 

i=i \ 



(1) 



where N is the system size, pi is the momentum of the 
ith particle, qi is its displacement from the equilibrium 
position, U is the internal potential, and V is an external 
potential. In the following we shall always use m — I, 
without loss of generality. Real isolated systems must 
have V — 0, since any external potential would itself in- 
volve a "scaffolding" which is not rigid and would thus 
also contribute to heat conduction. Thus, for real sys- 
tems the support of the external potential must be in- 
cluded in the description of the system, and V has to be 
replaced by a contribution to the internal potential U. 
Nevertheless, we shall in the following keep the ansatz 
(1), understanding that we are calculating only part of 
the complete heat conduction when V 0. In the follow- 
ing we shall also assume that the system has no frozen 
disorder. 

It is well known that heat conduction is infinite, if all 
potentials are harmonic 2] or if the system is integrable. 
In this case phonons resp. solitons are not scattered. 
Thus they propagate ballistically, given a constant heat 
flux J (independent of N) when a finite constant tem- 
perature difference AT is applied to the two ends of 
a chain of length N. Thus formally, the conductivity 
K — JN/AT is proportional to N. For non-integrable 
models with an acoustic phonon branch, i.e. without 
an external potential, one expects ballistics transport in 
the infinite wave length limit (since phonon scattering in 
general decreases with energy), which leads to a power 
behavior n ^ with < a < 1 (sj- This is the case, 
e.g., for the Fermi-Pasta-Ulam (FPU) model fij and for 
the diatomic Toda lattice . 

An exception to this seems to happen for coupled ro- 
tators nia, 



which seem to have finite k @ . This is explained by 
the fact that single highly excited rotators essentially de- 
couple from neighboring rotators, acting thus as barriers 
for the propagation of any phonons, even soft ones. 

Soft acoustic phonons are essentially Goldstone modes 
due to Galilei invariance. Finite heat conduction (in 
d = 1) is therefore expected when V ^ 0, since then 
translation invariance is broken and an acoustic branch 
does not exist in the phonon spectrum. This is e.g. the 
case for the Frenkel-Kontorova (FK) model 
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(notice that this is the commensurate FK model, where 
the harmonic part of the potential leads to the same par- 
ticle distance in the ground state as the cosine potential; 
we shall only discuss this case in the following). In this 
model, all phonons have a finite mean free path, bounded 
from above by a finite constant which is independent of 
the wave number k, but which diverges for T and for 
T oo. The latter follows from the fact that V given by 
Eq.(3) effectively becomes negligible when T — > oo, and 
effectively becomes a sum of harmonic potentials when 
T ^ . Thus there is no ballistic transport, and no 
obvious mechanism which could lead to an infinite con- 
ductivity for any finite T, while one expects k to diverge 
when T— s-OorT— >oo 11]. 

It was thus very surprising when Savin and Gendelman 
[T^ claimed to have clear evidence for phase transitions 
in the FK model, at which the conductivity changed from 
finite to infinite. They claimed that k is finite only in an 
interval Td < T < Tc2, with Td and Tc2 dependent on 
e, while k = oo outside this interval. 

Indeed, the same authors had also claimed that there 
is a phase transition in the rotator model 0, ^3 ■ There, 
the density of highly excited rotators should of course go 
to zero for T — > 0. Thus one expects that soft phonons 
exist in this limit, and k oo for T — > 0. Instead of this, 
it was claimed in |^ Il0| that k = oo in an entire interval 
0<T <Tc. 

In the present paper we want to test these claims by 
performing simulations on larger lattices and with higher 
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FIG. 1: Temperature profiles for tlie FK model with e = 1 and 
Thigh = 0.55, Tiow ~ 0.45. The lengths of the central parts 
of the chains are 64, 128, 256, 512, 1024, and 2048. While 
the gradients in the central region are roughly constant, most 
of the temperature variation happens for short chains in the 
thermostated boundaries. 



precision than in 0, ^| . 

In order to mimic the simulations of ^ Hol IT^ as close 
as possible, we also used Langevin thermostats (we do not 
agree with these authors that Nose-Hoover thermostats 
would be unsuitable, but we just don't want any discus- 
sion about this point). More precisely, we simulated a 
chain of A^o + + oscillators. The central TV os- 
cillators follow their Hamiltonian equations of motion, 
while the outer 2No ones satisfy gr„ = — ^ — jQn + 
and with (p — U + V being the total potential, ^„ being 
white Gaussian noises, {^n{t)^kit')) = 2jTnSnkS{t — t') 
and with T„ = Thigh for — A^o <n<Q and T„ = Tiow for 
N <n< N + Nq. We used No = 40 and 7 = O.I, as in 
duo, El- The temperature difference AT = Thigh -Tow 
was chosen between 10% and 20%. 

For the integration we used a simple leap frog 'l^- 
On the one hand this is symplectic and thus more suited 
for the central region than, say, a Runge-Kutta integra- 
tor. On the other hand it should be more robust than 
higher order symplectic integrators in the boundary re- 
gions which are not Hamiltonian. Step size was 0.05, and 
total integration times were typically ~ 10^ — 10^ units 
(i.e. 10® — 10^ steps), with some runs going up to 5 x 10* 
units. We checked that this was sufficient to reach a 
steady state and that the time-averaged heat flux J was 
independent of the site. 

We verified also that the temperature profiles were 
roughly linear in the central region < n < TV, but 
we could not verify the absence of temperature jumps at 
its boundaries claimed in ^ llolIT^ . More precisely, such 
jumps were absent only for very large lattices and small 
conductivities, i.e. if the heat flux was small. Other- 
wise, for small lattices and/or large conductivities, there 
were very large jumps, mostly located in the boundary 
regions -40 < n < and < n < -I- 40 (see Fig. I). 
Thus while the profile in the central region was essentially 
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FIG. 2: Heat conductivity versus system size for the model. 
Each curve corresponds to a constant value of /3, while T = 2 
is held fixed. Statistical and integration errors are less than 
the symbols. The continuous lines are only drawn for guiding 
the eye. 



linear (in contrast to simulations with Nose-Hoover ther- 
mostats coupled to single particles, as e.g. in 0,0], it 
would be very wrong to estimate the conductivity sim- 
ply by dividing the flux by the nominal imposed AT. It 
seems that this was done in several cases in ^2 ^ 

which explains some - but not all - of the differences 
between their results and those of the present paper. In 
other cases the authors of 0,0,0 must have taken into 
account boundary jumps, otherwise their results would 
disagree much more with ours than they actually do. Un- 
less otherwise said, we will estimate k, by dividing the flux 
by the temperature drop over the inner half of the central 
region. 

Notice that a similar behavior was found also for the 
rotator model and for the discrete (f)^ model. The latter 
is given by the Hamiltonian )14j 
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pi + (ft+i - lif + c^Qt + ^Pqf 



(4) 



Conductivities for the (p'^ model with T — 2 and (3 = 1 — a 
are shown in Fig. 2 for various values of /3 € [0, 1]. The 
0^ model is a prototype model with finite conductivities. 
Indeed we see that all measured values of k, are not only 
finite but are independent of A^. This would not have 
been the case if we had not taken the temperature jumps 
into account and would have used the nominal value of 
AT when estimating n from J/ AT. The values of k 
shown in Fig. 2 are in very good agreement with those of 

Frenkel-Kontorova model: In all simulations, A^ 
ranged from 32 to 2048. We made simulations for e — 
1.0, 3.0 and 10.0. According to 0| , Td is only weakly de- 
pendent on e: = 2.6, 2.3, and 2.0 for e = 1.0, 3.0, 10.0. 
On the other hand, Tc2 should strongly increase with e, 
Tc2 ~ 3.3, 15, and 150 for the above three values. 

Conductivities for e = 1.0 are shown in Fig. 3. Ob- 
viously, they are finite for T = 2 and T = 8, showing 
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FIG. 3: Heat conductivity in the FK model versus system size, 
for e = 1.0. Each curve corresponds to a constant average 
temperature. 
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FIG. 4: Heat conductivity versus system size for e = 3.0. 
Statistical and integration errors are less than the symbols. 



that the temperature range with finite conductivities is 
underestimated in 12]. For T = 0.5 and T = 32 we see 
a slow increase of k with TV, over a wide range of the 
latter, but it seems to stop for the very largest lattices 
{N > 1000). Finally, for T = 0.125 and T = 128 there 
is a slow increase for all N > 300. The latter could be 
taken as an indication that k diverges for these temper- 
atures, but we think that this would be wrong. On the 
one hand, the increase with N is very slow, much slower 
than in Fig. 4 of We would get a similarly fast in- 
crease as in if we would use the nominal temperature 
difference, i.e. if we would disregard the jumps seen in 
Fig. 1. On the other hand, the data for 0.5 < T < 32 
show us that the saturation of k happens at larger and 
larger lattice sizes as we go away from the central energy 
region, just as we have expected. Thus we must expect 
in any case an increase of k for all reachable lattice sizes, 
and observing it does not give any relevant information. 

Analogous results for e = 3.0 are given in Fig. 4. There 
we show data for the temperature range [0.75,75]. This 
is again much larger than the range where convergent 
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FIG. 5: Asymptotic (for N — » oo) heat conductivity versus 
temperature for e — 3.0. The x axis is the temperature, the 
y axis is conductivity Kccntor- Actually, the plotted values of 
K are those measured for N — 2048, but they seem to be 
independent of lattice size for N > 1000. 
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FIG. 6: Heat conductivity versus system size for e = 10.0. 
Again, statistical and integration errors are less than the sym- 
bols. 



conductivities were found in 5l2|. This time all curves 
become horizontal for large i.e. the conductivity is 
finite in the entire range. It of course depends strongly 
on T, see Fig. 5. It diverges both for T ^ and for 
T — > oo, since the problem effectively becomes harmonic 
in both limits. 

Finally, our last simulations for the FK model, for e = 
10, are summarized in Fig. 6. There we only show results 
for low temperatures, 0.5 < T < 2.75. Except for the last 
temperature, they are all in the regime where the authors 
of (12i] have found divergent n. In contrast, all our curves 
are either horizontal for all N or become horizontal for 
large N, suggesting that k is finite for all finite T. 

Rotator model: For the rotator model, we simu- 
lated larger systems, with ranging from 32 to 8192. 
Conductivities are plotted in Fig. 7 against A^ for T = 
0.6,0.45,0.3, and 0.2 (from bottom to top). According 
to the phase transition from a high-T phase with 
finite conduction to a low-T phase with infinite conduc- 
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FIG. 7; Heat conductivity versus system size for the rotator 
model. Again, statistical and integration errors are less than 
the symbols. 

tion occurs at some Tc between 0.2 and 0.3. According 
to that, the lowest three curves in Fig. 7 should become 
flat for TV — > cxD, while the uppermost should continue to 
grow. This is not what is found, although our values for 



T = 0.2 and 0.3 agree numerically quite well with those 
of But while the curve for T = 0.2 increases with 
the same average slope as in Fig. 2 of , it is definitely 
S-shaped and stops to rise for the largest values of N . 

Conclusion: In this paper we studied the size depen- 
dence of the effective finite size conductivity of nonlin- 
ear ID lattices, as a function of temperature. We used 
straightforward but high statistics simulations to show 
that there are no indications of the phase transitions sug- 
gested in 1^ ITol I12I I on the basis of similar simulations. 
For the FK model, this is in agreement with expectations, 
since phonons should have a finite free path in this model 
for all finite temperatures. For the rotator model it is less 
obvious. It suggests that the blocking of the propagation 
of soft phonons by localized excitations is effective at 
all finite temperatures. It becomes of course less and less 
important as T ^ 0, since the density of such excitation 
decreases exponentially with 1/T. But it is present at all 
finite T, and it becomes dominant for N 00. 

P. G. wants to thank Roberto Livi and Antonio Politi 
for very helpful discussions. 
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